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_ ABSTRACT 

^D I By fitting the classical Rotating Vector Model (RVM) to high quality polarization data for selected 

C^ ■ radio pulsars, we find the inclination of the magnetic axis to the spin axis, a, as well as the minimum angle 

between the line of sight and the magnetic axis, f3, for ten objects. We give a full treatment of statistical 

^ |. errors in the fitting process. We also present a dictionary and conversion table of various investigators' 

QJ • geometric definitions to facilitate future comparisons. We compare our results with other RVM fits and 

C^ I with empirical / geometrical (E/G) approaches, and we examine the strengths and weaknesses of RVM 

QQ . fits and E/G investigations for the determination of pulsar emission beam geometry. 

Our fits to B0950+08 show that it is an orthogonal rotator with the main and interpulse radiation 
emitted from opposite magnetic poles, whereas earlier RVM fits indicated that it is an almost-aligned, 
single-magnetic pole emitter. We demonstrate that low-level emission across a wide longitude range, 
K* ' when properly weighted in the RVM fit, conclusively favors the former scenario. B0823+26 is also 

^«0 , an orthogonal rotator. We find that B1929+10 emits into its wide observed range of longitudes from 

j^ ■ portions of a single cone that is almost aligned with the spin axis. This result agrees with other RVM 

fits but conflicts with the E/G findings of Rankin & Rathnasree (1997). 

We determine that convergent RVM solutions can be found only for a minority of pulsars: generally 

^^ ' those having emission over a relatively wide longitude range, and especially those pulsars having inter- 

f— ^ I pulse emission. In pulsar B0823+26, our preferred fit to data atall longitudes yields a solution differing 

by several a from a fit to the main pulse / postcursor combination alone. For pulsar B0950+08, separate 

fits to the main pulse region, the interpulse region, and our preferred fit to almost all longitudes, converge 

to results differing by several times the formal uncertainties. These results indicate that RVM fits are 

Q ■ easily perturbed by systematic effects in polarized position angles, and that the formal uncertainties 

H I significantly underestimate the actual errors. 

iG ' Subject headings: pulsars: geometry - pulsars: magnetic inclination - pulsars: interpulses - pulsars: 

rotating vector model 
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1. INTRODUCTION derive and apply our own technique for determining pulsar 

5-h' t-,, •,• f-.i , ri- 1 •• orientation angles, and we compare our method with oth- 

rri • Determmation 01 the geometry of radio pulsar emission ,. , , '=., ' ...,,,■, , , r 

^^ , ■ i- 1 i J i J- • • \ ■ rni. ers m order to illuminate the strengths and weaknesses 01 

■IS essential to understanding emission mechanisms. Ihe , ,, ° 

orientation of any pulsar reduces basically to two angles: 
a, the angle between the spin axis and the observable mag- 
netic axis, and f3, the minimum angle between the mag- IT. The Rotating Vector Model (RVM) 
netic axis and the observer's line of sight as the beam One way to obtain a and /3 is to examine the position an- 
sweeps past the observer (see Fig. 1). Finding values for gj^ ^f ^^^ linearly polarized emission from the pulsar, and 
these angles can lead to the determination of the intrinsic -^-^^^^^ ^^-^ -^ ^^^ method that we will use below on our 
beam width and other geometrical properties of the pulsar ^^^ ^^^^^^ The S-shaped curve visible in plots of the polar- 
emission. For example, knowing a for a range of pulsars j^ation position angle vs. longitude (see, for instance, the 
gives us clues about the origin of their magnetic fields and ^^^^^^ p^^^gj of pig^ 3) jg explained by a model proposed 
how the pulsars are evolving (Candy & Blair 1986; Biggs ^^^^^■^y ^^^^^ ^^^ discovery of pulsars (Radhakrishnan & 
1992; McKmnon 1993), while knowledge of /3 leads to ^ooke 1969) in which the position angle follows the rota- 
mformation on beam properties (Biggs 1990; McKmnon ^^^^ ^f ^^^ magnetic field lines at the sub-Earth point on 
1993). In addition, theories about the structure of pul- ^he pulsar. Called the Rotating Vector Model (which we 
sar beams make predictions for the orientation of partic- j^note hereafter as RVM), the model gives the polariza- 
ular pulsars based on empirical and geometrical equations ^-^^^ position angle, V', as a function of pulsar longitude, 
withm each given theory, so independent methods can be ^^ ^^^^^ ^^^ p^jg^^ rotation equals 360° of longitude: 
used to verify those predictions and hopefully help dis- 
tinguish between the models. In this paper, we carefully , / ,/ ;'\ — sinasin(0 — (j)o) 

ian(ip — ipQ) — — — ——. -— -—■. [1) 

smCcosa — cosCsinacos((/) — <po) 

^Also Department of Astronomy and Astrophysics, University of Chicago, 5640 S. Ellis Street, Chicago, IL, 60637. 
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(We use "0 for the observed polarization position angle, 
which will be shown below to be different from ip' , below.) 
The offset angles -00 a-nd 0o (constant for each pulsar) 
give the polarization position angle and longitude, respec- 
tively, of the symmetry point and maximum gradient of the 
position angle curve (when the pulsar beam is pointed clos- 
est to the observer), and also include arbitrary constant 
offsets associated with observing parameters. The quan- 
tity a is the angle between the positive spin axis (which 
points in the direction of the angular velocity vector il) 
and the observable magnetic pole, while Q is the angle be- 
tween the positive spin axis il and the pulsar-observer line 
of sight. The sign and magnitude of /3, the impact param- 
eter of the line of sight with respect to the magnetic axis, 
are defined by 

/3 = C-a. (2) 

In the context of these definitions, (3 is positive whenever 
^ > a so that the line of sight vector is farther from the 
(positive) spin vector +Vl than is the observable magnetic 
axis; and j3 is negative whenever C < a so that the line of 
sight vector is closer to the (positive) spin vector -\-Q. than 
is the observable magnetic axis. In our RVM fits that fol- 
low, we hold to these definitions regardless of whether a is 
greater or less than 7r/2, so that we remain true to the con- 
ventions of Eq. 1. It is important to note that a positive j3 
corresponds to an "outer" (i.e., equatorward^) hne of sight 
only if a < 7r/2; whereas for a > 7r/2 the equatorward line 
of sight has a negative (3. Other investigators use differ- 
ent conventions for the sign of /3, as we will discuss below. 
Furthermore, as pointed out by Damour & Taylor (1992) 
and Arzoumanian et al. (1996), Eq. 1, used by essentially 
all researchers who have fit the RVM model to data, was 
derived with the convention that the polarization position 
angle, tp' , increases clockwise on the sky. This is contrary 
to the usual astronomical convention that measured polar- 
ization position angle "0 increases counterclockwise on the 
sky. Since most previous analyses have fitted the RVM via 
Eq. 1 and its clockwise ip' to data having the observers' 
convention of counterclockwise ip without correction, we 
have modified many earlier investigators' results in order 
to be consistent with the assumptions of the RVM model. 
In what follows, we refer to this issue as the "'0 conven- 
tion problem." An RVM fit with the "0 convention problem 
must be corrected by transforming the published values, 
which we refer to as aorig and (3orig, to values a and /3 
that adhere to the clockwise-'0' RVM rotation convention 
of Eq. 1, which we also use in our fitting process: 



magnetic pole. The slope of this sweep contains impor- 
tant geometrical information that is used either explicitly 
or implicitly by essentially all investigators: 



a = 180° - a 



ortg 



(3) 
(4) 



We will consider this conversion in more depth when we 
compare our results to those from other workers. At that 
point we will also show that definitions of some other quan- 
tities in some earlier works must also be modified for con- 
sistency with Eq. 1 and its associated definitions (see Ta- 
ble 1 for conversion relationships). 

The ip convention problem manifests itself in another 
manner as well. One of the most characteristic properties 
of the RVM is the steep swing of polarized position an- 
gle as the line of sight makes its closest approach to the 

^ "Equatorward" indicates that the Une of sight is opposite the spin 
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Note that since the observers '-defined ip = —ip', a minus 
sign must be inserted in the above equation if ip' is re- 
placed by -0. Note also that Eq. 5 fixes the sign of /3 since 
sin a > 0. However, as emphasized in the discussion fol- 
lowing Eq. 2, the sign of (3 does not by itself select inner 
versus outer line of sight trajectories. See Table 2 for a 
summary. 

1.1.1. The Second Magnetic Pole 

The second magnetic pole is relevant in some cases of 
interpulse emission. The opposite magnetic pole's colati- 
tude is a2 = tt — a, and its impact parameter with respect 
to the line of sight is /32 = C ^ 0:2- It is important to note 
that the RVM itself does not distinguish one- from two- 
pole emission, as the model provides the position angle of 
the particular magnetic field line that happens to be at 
the sub-Earth point at any instant, which is a function of 
magnetic dipole geometry alone. 

1.2. Earlier RVM Fits 

Many researchers have attempted to fit the RVM func- 
tion to polarization data (the most comprehensive of which 
include Narayan & Vivekanand (1982); Blaskiewicz et. al. 
(1991); von Hoensbroech & Xilouris (1997a, b)) in order to 
determine a and /?. In comparing those results to our data, 
however, we have found that there has been a wide range 
of definitions of geometrical beam angles among the dif- 
ferent authors. We now discuss their procedures in detail. 
In order to compare the different results, we also present 
for the first time a dictionary and procedure for converting 
among the various investigators' definitions (see Table 1). 
The earlier results, rationalized by the rules of Table 1, 
are shown in Table 3 along with our results. 

1.2.1. Narayan & Vivekanand (1982) (NV82) 

Narayan & Vivekanand (1982) fitted the RVM model 
to single-pulse polarization histograms, thereby eliminat- 
ing problems caused by emission of orthogonal polarization 
modes (see below and Backer & Rankin (1980)). They un- 
weighted those longitudes where the position angle could 
be determined in less than 15 percent of the pulses, and 
uniformly weighted the rest. They emphasized the diffi- 
culty of distinguishing outer (i.e., equatorward) and inner 
line-of-sight trajectories from RVM fits when the pulsar's 
emission occupies only a few percent of the pulse period, 
as is usually the case. Despite this difficulty, their fits did 
tend to favor one trajectory over the other in most cases. 
Their a was measured with respect to the nearest spin 
pole so that it is never greater than 7r/2 in their full RVM 
fits, and they assigned outer (equatorward) lines of sight 
to positive j3 in all cases. These definitions, eminently 
defensible on physical grounds, are nonetheless at odds 
with those of Eq. 1. Substantial gymnastics are needed 
to transform from one to the other. Uncertainties were 

pole lying nearest to the observable magnetic pole. 
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given for only some of the results, and it was not men- 
tioned whether the strong covariance of a and (3 (which 
we find in our own fits, and elaborate on below) is reflected 
in these uncertainties. 

1.2.2. Blaskiewicz et al (1991) (BCW91) 

These authors fitted the RVM to a subset of the Wcis- 
berg et al. (1999) Arecibo 1418 MHz pulsar polarization 
data and also searched for manifestations of special - rcl- 
ativistic effects. Publishing while the experiment was still 
in its data-acquisition phase, they necessarily used shorter 
total integrations than the data analyzed and displayed 
here. They chose to exclude data with L below 5 or 10 
times its off-pulse RMS as well as regions appearing to 
contain orthogonal modes (see Backer & Rankin (1980)). 
Also, they used uniform weighting for the data that sur- 
vived to be fitted. In the end, out of 36 fits on 23 pulsars at 
0.43 and 1.4 GHz, only 17 had uncertainties in a that were 
< 50°. Fractional uncertainties in (3 were often also very 
large. Again, it was not mentioned whether the covariance 
of a and (3 is reflected in these errors. Their solutions are 
reported here after correcting for the if) convention prob- 
lem. We will see below that many of our results have 
smaller uncertainties than BCW91, presumably because 
of our larger quantity of data and consequently improved 
signal-to-noise ratios, our use of non-Gaussian statistics 
and non-uniform weighting in the fitting procedure. 

1.2.3. von Hoensbroech & Xilouns (1997a,b) (HX97a,b) 

More recently, von Hoensbroech & Xilouris (1997a, b) 
also attempted RVM fits to polarization data at 1.4, 1.7, 
and 4.85 GHz. They used the Simplex algorithm to ap- 
proach the x^ minimum and then zeroed in with the 
Levenberg-Marquart fitting algorithm, much as we do be- 
low. Their quoted uncertainties are also generally rather 
large, which they attribute partially to the high correla- 
tion of a and /3. Goodness of fit was assessed with the x^ 
test, with a study of the departure from Gaussian statis- 
tics of the post-fit residuals, and through evaluation of the 
symmetry of post-fit residuals around the best fit. There 
is no mention of a threshold noise cutoff nor of differen- 
tial weighting of the points. In Tables 1 and 3, wc correct 
their results for the ^p- convention problem, and for cases 
having their anx < 0. 

1.3. Previous Empirical / Geometrical (E/G) Solutions 
for a and (3 

Combinations of empirical and geometrical approaches 
have also been used to help determine beam alignments. In 
these cases, a full fit to the polarized position angle curve 
is not performed. Instead, certain empirical relationships 
are derived and combined with beam geometry equations. 
The resulting sets of equations are then evaluated as a 
function of observed pulsar properties. By their nature, 
these techniques do not yield estimates of uncertainties. 

1.3.1. Lyne & Manchester (1988) (LM88) 

Lyne & Manchester (1988) determined a and /3 for pul- 
sars which they believe exhibit emission from both sides 
of a circular cone, having a full longitude width at 10% 



maximum of 2A(/i. It can be shown geometrically that the 
emission cone intrinsic angular radius, p, is given by 

sin ( — ) = sin ( ) sinQ:sin(Q! + /3) + sin ( — ). (6) 

LM88 calculated hypothetical perpendicular rotator pa- 
rameters "/Jgo" and "/Ogo" for each such pulsar, where 

/390 = P\a=90° ^'^'^ ^90 = PL=90° • F''"™ ^q. 5, 

sin/390 = [%\„ra.]-^\ and Eq. 6 yields sin\l^) = 

sin^(^)cos/39o -I- sin^(^). A plot of pgo as a function 
of pulsar period, P, reveals a lower limit to the scatter of 
data points. LM88 argue that the lower limit represents 
pulsars that truly have a = 90°, in which case the intrinsic 
beam radius p = pgo, whereas the points scattered above 
the limit represent the case a ^ 90° so that the true beam 
width is p < pgo- From the lower limit, they derive the 
important result: 

p ^ 6°5P-^/^. (7) 

Calculating p from this relation, and using the other mea- 
sured quantities dtp' / d<j)\max and 2A(j), it is then possible to 
use (the unapproximated) Equations 5 and 6 to calculate 
a and /?. 

These calculations rest on the assumption that the pul- 
sar beam is circular (which is backed up by the favorable 
comparison of their calculations for the range of the po- 
sition angle swing (2Ai/') and their measurements of that 
swing), as well as their idea that certain pulsars have "par- 
tial cone" emission for which the above relationship cannot 
be used. Concerning the first assumption, there has been 
much controversy over the shape of pulsar beams, with 
some investigators advocating emission elongated in the 
latitude direction (Jones 1980; Narayan & Vivekanand 
1983a), others finding extension in the longitude direction 
(Biggs 1990), (McKinnon 1993), and yet others agreeing 
with LM88 that the beam is essentially circular (Bjorns- 
son 1998). It should be noted that later work by this 
group (Manchester, Han, & Qiao 1998; Tauris & Manch- 
ester 1998) suggests a different functional form for Eq. 
7, which if applied to the pulsars considered here, would 
yield somewhat different values for a and (3. Also, it has 
been argued by Rankin (1990) that some of the types of 
pulsar emissions under consideration here could be what 
she calls "core" beams, for which she finds the RVM, and 
hence Eq. 5, problematic. 

The results from Lyne & Manchester (1988) were very 
helpful to us as initial guesses for our own fits. However, 
these investigators tabulated only \P\ and confined a, by 
definition, to 0° < a < 90°. We defined the sign of their 
f3 from the sign of the position angle sweep (c/. Eq. 5), 
but for a we had to choose between ulm and tt — a^M 
solely on the basis of their consistency with our final re- 
sults. Note that this represents the choice between inner 
and outer lines of sight trajectories, so it is not trivial (see 
Section 1.1 and Table 2 for further discussion of this issue). 

It should be noted that LM88 also provided, in addition 
to their E/G analyses, RVM fits to several pulsars having 
interpulse emission. We will discuss these results in the 
sections on those individual interpulsars. 

1.3.2. Rankin (1990, 1993a,b) (R90, R93a,b) 

In a series of papers, Rankin (1983, 1990, 1993a, b) laid 
out a comprehensive pulsar classification scheme that dis- 
tinguishes two basic types of emission: core pencil beams 
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and hollow cones. R90 discovered a simple relationship 
between intrinsic core pulse width and period by studying 
pulsars thought to have opposite-pole interpulses. Sonic 
simple geometry then permits the calculation of a for any 
pulsar with a core component. In this method, Wcore, the 
longitude FWHM of an interpulsar with a core component, 
was measured and interpolated to 1 GHz. In a remarkably 
good fit, six intcrpulsars all follow the relation 



Wr.. 



= 2? 45 P 



-1/2 



(8) 



Presumably then for these pulsars, which must have 
a ~ 90° in order that both pulse and interpulse emis- 
sion from opposite poles be observable, Wcore measures 
the intrinsic width of the core beam. All other core pul- 
sars have larger observed Wcore for a given period, which 
is consistent with the idea of a beam of the same intrinsic 
width having a larger longitude extent Wcore if directed 
away from the equator (i.e., a ^ 90deg): 



Wc, 



2?45 p-^^^/sina. 



(9) 



R90 then took measured values of Wcore and P to solve 
for a via Eq. 9. Strictly speaking, this procedure yields 
only sin a, so that it is not possible to distinguish a from 
TT — a. It should also be noted that this equation assumes 
that the impact parameter |/3| is small, and ignoring /3 is 
defended on the grounds that its effect on core compo- 
nents is weak because the angular intensity may be ap- 
proximated by a bivariate Gaussian, whose FWHM is in- 
dependent of (3. It is also assumed that the core emission 
beam completely fills the open field lines, so that the core 
width is directly related to the angular width of these field 
lines. R90 did not attempt to determine /3 from these core 
components, stating that "the polarization-angle behavior 
of core components seems to provide no reliable informa- 
tion about the impact angle (3." 

In a later pair of papers, Rankin (1993a, b) did use po- 
larization position angles to determine /3 for pulsars having 
conal emission. The position angle swing relationship, Eq. 
5, was used, so both dip' / d(j)\max and a must be deter- 
mined. The slope was determined from the observed posi- 
tion angle curve. For pulsars having both core and conal 
emission, a was calculated as above, while other tech- 
niques were used in the absence of core emission. R93a,b 
attempted to distinguish outer and inner line of sight tra- 
jectories by noting whether the position angle curve did or 
did not flatten in the wings (c/. Narayan & Vivckanand 
(1982)), but this was not based on rigorous fits. Her sign 
convention then followed NV82 (positive /3 for outer line 
of sight trajectories), although we found that we occasion- 
ally had to flip her published signs for internal consistency 
(see Table 3). 

Interestingly, R93a,b determined that cones, like cores, 
also have intrinsic widths that depend only on period. 
Two angular radii are indicated, suggesting the existence 
of both an inner and an outer cone (not to be confused 
with the inner and outer line of sight trajectories discussed 
above) . 

1.4. Summary 

The various efforts at RVM fits involve subtle assump- 
tions regarding choice and weighting of data. Both of the 



empirical / geometrical methods also involve varied as- 
sumptions, many of which reach to the heart of theories 
of pulsar emission. Interestingly, although the two E/G 
techniques are based on different assumptions and dif- 
ferent definitions, they generally yield similar results for 
a < 40°(Miller & Hamilton (1993)). 

We have fitted our own time-averaged data to the RVM 
to assess the efficacy of both empirical-geometrical and 
RVM approaches to determining the pulsar and beam ori- 
entation. Work on fitting single-pulse data is ongoing, and 
will be addressed in a separate paper (Weisberg & Everett 
2001). 

2. OBSERVATIONS 

The average pulse data we analyze here were taken at 
Arecibo Observatory at 1418 MHz from 1989 through 1993 
(Weisberg et al. 1999, (W99)). For these data sets, we 
attempted to fit the RVM to every pulsar whose position 
angle appeared to follow the RVM model. We were able 
to find rigorously convergent fits on only a small fraction 
of the pulsars in the database, for reasons that wc will 
discuss below. 

3. OUR FITS 

In completing our fits to RVM polarization data, we 
found that wc had to be especially careful in deciding on 
statistical vs. uniform weighting, calculating errors in our 
polarization data, dealing with the possibility of orthogo- 
nal modes, fitting covariant parameters in the RVM, and 
calculating the errors on those covariant parameters. Each 
of these concerns is addressed in the subsections below. 

3.1. Statistical Vs. Uniform Weighting 

One of the most important limitations of the RVM fit- 
ting technique is that a large longitude range of polariza- 
tion position angles is necessary in order to find trustwor- 
thy values for a and /3 (Narayan & Vivckanand 1982; 
Miller & Hamilton 1993). The narrowness of most pul- 
sar beams unfortunately renders this goal difficult in most 
cases. Thus, many of the best RVM fits are to pulsars with 
interpulses, where polarization information is available at 
a very large range of longitudes (e.g., Lyne & Manchester 
(1988); Phillips (1990)). This difliculty is aggravated by 
the fact that most investigators have uniformly weighted 
all data used in an RVM fit, eliminating any that falls 
below some sensitivity threshold or outside of some longi- 
tude range, and therefore reducing further the amount of 
polarization information remaining to be analyzed. In or- 
der to use as much data as possible and therefore decrease 
the uncertainty in our fits as much as possible, we im- 
plemented a statistical weighting scheme. This procedure 
enables us to sample a larger range of longitudes and to 
use xt to asses the quality of the fit. However, to use sta- 
tistical weighting to its fullest advantage, one must be very 
careful in calculating uncertainties for each data point. It 
is also important to bear in mind that systematic errors 
cannot be accounted for in this scheme. 

3.2. Error Calculations for Polarization Data 

For the Stokes parameters /, Q, U, and V, each asso- 
ciated standard deviation obeys Gaussian statistics, with 
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o'l = aQ — au — ay- However, polarized quantities de- 
rived from the Stokes parameters such as linearly polarized 
power L and position angle tp are not normally distributed, 
and L itself suffers a bias. In order to correctly use polar- 
ization data, we must take these issues into account. We 
show first how we removed the bias in L, and then we dis- 
cuss the non-Gaussian errors in tp. (Errors in tp and in tp' 
are identical, so for simplicity wc will refer only to a^ in 
what follows.) 

3.2.1. Removing the Bias from the Measured Linear 
Polarization Lmeas 

To compute Lmeas, we used the standard definition: 

Lmeas = VQ^TIP . (10) 

The Ljrieas Calculated in this way is biased (Naghizadeh- 
Khouei & Clarke 1993), however, because it is defined as 
a positive definite quantity and the simple equation above 
overestimates the true polarization in the presence of noise. 
This can be seen by noting that if the true values of Q and 
U in the above equaton for L were both zero, noise on both 
of those values would lead to a positive value for Lmeas- 
Thus, the error in the Stokes parameters always increases 
the value of L, and it is this non-zero bias that we wish to 
remove from our value of Lmeas ■ 

Simmons & Stewart (1985) rigorously compared the 
merits of various schemes for debiasing L, and showed that 
several of the estimators lead to negligibly different out- 
comes, especially at intermediate signal-to-noise ratios. 
We chose the estimator of Wardle & Kronberg (1974) 
for all but the lowest signal strengths, where we set the 
unbiased linear polarization to zero following Simmons & 
Stewart (1985). We first calculated the off-pulse standard 
deviation of the total power Stokes parameter /, aj. The 
measured quantity Lmeas /o'l, which can be considered a 
kind of "signal-to-noise," then serves as the correction pa- 
rameter in the expression for finding the unbiased linear 
polarization Ltme'- 



v^ 



'-Y - 1 fj/ if -^ > 1.57 



otherwise. 



3.2.2. Error Estimation for Polarized Position Angle tp 

At high signal-to-noise (S/N) ratios, the uncertainty in 
the position angle is easily calculated: 



., . |i . 28°65 ^ 



(11) 



However, because of the above-discussed bias to the linear 
polarization, repeated measurements of the position angle 
of polarization at low- to intermediate-S'/TV are not nor- 
mally distributed, and special formulae must be used to to 
determine a^. The probability distribution of the position 
angle (tp) around the true value (tptrue) can be calculated 
(Naghizadeh-Khouci & Clarke 1993): 



G{tp; tptrue, Po[ 






ryoe''"[l-herf(77o)]}e 



-i^) 



where rjo — -^ cos 2{tp — tpi 



Pn = 



(12) 
(note then 



that L must be debiased as shown above before Pq can be 
calculated), and "erf" is the Gaussian error function. 



Following Naghizadeh-Khouei & Clarke (1993), we de- 
fine "Icr^" via the numerical integral of the above function 
(with tptrue = 0) by adjusting the bounds of integration to 
include 68.26% of the distribution, as below: 



G{tP,Po)dtP=^ 68.26% 



(13) 
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To find the error on any given data point, we first built 
a table of the results of this numerical integration for Pq 
in the range of 0.0 to 9.99 in steps of 0.01. Then, to find 
the "Icr^" error bars for a data point having Pq < 10, 
we linearly interpolated on this table. At higher signal 
strengths, Eq. 11 was used. 

We used the values of a^ from this procedure as the 
uncertainties for each longitude bin in our fits (and thence 
as the weighting factors). 

3.3. Orthogonal Modes 

Pulsar radiation is frequently found to be emitted in two 
orthogonally linearly polarized modes at the same longi- 
tude. Variations in strength between the two modes can 
severely corrupt a classical RVM curve. To deal with this 
possibility, we used the following procedures: In those rare 
cases where a 90° position-angle discontinuity suggested 
that emission was switching from predominantly one mode 
to predominantly the other, we inserted 90° discontinu- 
ities into our model as well. In the more common cases 
of RVM disruption due to a slow change in the balance 
of two orthogonally emitted modes, we unweighted the af- 
fected longitudes in our fits. We identified these possible 
orthogonal modes by referring to single pulse studies when 
possible (Backer & Rankin 1980; Stinebring et al. 1984). 

3.4. Fitting Covariant Parameters 

We then fit the RVM model to tp{(p) (and a^{(p)) using 
mainly the Levenberg-Marquardt algorithm (Press et al. 
1992) due to its ability to work efficiently with non-linear 
functions and to yield good estimates of the errors on the 
parameter values it finds. 

During the development of the RVM fitting engine, an 
RVM simulator was written to ensure that the fitting code 
was working properly. In order to generate realistic simu- 
lated data, we first created a classic RVM position angle 
curve, and then added the non-Gaussian noise and biased 
L{(t)) [see above for discussion of the noise and biasing]. 
Our noisy, biased simulated tp{(p) and L{(p) were then fed 
into our standard fitting program. Multiple runs of the 
simulator and fitting program demonstrated that we could 
consistently find x^ k, 1±0.20, even with occasional (spu- 
rious) position angle data points far from the main pulse 
passing the S/N criteria. Comparisons of these advanced 
simulator runs with other editions of the simulator (with- 
out noise and debias code included) showed us that even 
considering the non-Gaussian distribution of position an- 
gle noise and biasing of off-pulse L noise, it should still be 
possible to use x^ as a reasonable goodness-of-fit indica- 
tor. 

For initial estimates for the parameters, we used our 
own best guesses as well as the results from Narayan & 
Vivekanand (1982); Lyne & Manchester (1988); Rankin 
(1990);Blaskiewiczet. al. (1991), and Rankin (1993a,b). 
As outlined in Sections 1.2 and 1.3, however, the results 
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given by other researchers for a and /3 (or whatever sym- 
bols they used for those angles) were not always consistent 
with the RVM equation, and had to be converted to a con- 
sistent definition (sec Tabic 1). 

Even with those initial guesses in hand, fitting the data 
is a very sensitive and difhcult process. The chief difficulty 
is that a and /3 are extremely covariant: without sufficient 
data in the "wings" of the polarization profile, the fitting 
engine is unable to lock into specific values for a and /3 and 
would find a -^ f80° and /3 ^ 0° (while still remaining 
consistent with the position angle sweep relationship, Eq. 
5). This problem is also aggravated when f/^o a-nd 0o are 
allowed to vary: the fitting engine can sometimes lock onto 
values for -00 and 0o and will then be unable to lower ^ 
without pulling (a,/?) toward (180°, 0°). These problems 
made fitting difhcult, but not impossible, at least in some 
cases. To help alleviate these problems, we added the Pow- 
ell fitting algorithm (Press et al. 1992) to refine our initial 
guesses before passing those results on to the Levenberg- 
Marquardt algorithm. We normally started by using the 
initial values from other researcher's results, and manually 
tweaking (^o and i/^o to fit our data. We then allowed the 
fitting engine to fit for a, /3, ^o and f/'o altogether. 

At other times, the fit would never truly converge to an 
answer, but would instead indefinitely search for a min- 
imum without finding one. The reason for the difhculty 
in fitting can be seen in a plot of the ^ topography of 
a particular pulsar, P0301+19, in Fig 2, which is repre- 
sentative of most. (See also von Hocnsbroech & Xilouris 
(1997a).) There is a deep chasm in y^ space representing 
all (a, /3) satisfying the observed position angle sweep via 

Eq. 5. Consequently the best-fitting solution (afcestj/Sfcest) 
lies somewhere in the bottom of the chasm. Unfortunately, 
the floor of the entire valley lies at approximately the same 
"elevation," resulting in a difhcult-to-find minimum and 
hence a poorly defined best fit. Note then that it is even 
difficult to distinguish something as basic as inner and 
outer trajectories. This is a graphically-oriented descrip- 
tion of the fact that a and /3 are highly covariant as a 
result of the similarity of all RVM curves satisfying a par- 
ticular value of Eq. 5, which is the limiting case of Eq. 1 
over the small longitude range occupied by most pulsars' 
emission. 

Another problem is that many pulsars seem to have 
RVM-style polarization curves, but upon closer inspec- 
tion, we discovered that the curve was of a slightly incor- 
rect shape, and again, the fit would not converge. Due to 
difficulties like these, out of 70 pulsars for which we had 
average pulse position angle data, only ten could be fit re- 
liably. The reliability of the fit was assessed by using other 
researchers' results as initial parameter estimates for our 
algorithm. In each of the ten successful cases, all of our 
fits converged to the same final results, which often were 
not the same as the other published solutions. The results 
for our fits are listed in Table 3, along with values from 
other researchers for those same pulsars, which have all 
been converted to consistent definitions for a and /3 via 
Table 1. 

3.5. Calculating Errors on the Covariant Parameters 

With the output from the fitting engine, we were able to 
calculate errors based on the full length of the region ly- 



ing at Ax^ = 1 above the minimum, which is imperative 
because a and /? are highly covariant (von Hoensbroech 
& Xilouris 1997a; Press et al. 1992). The estimated 
uncertainties reported in Table 3 reflect the covariant re- 
lationship between parameters. 

3.6. Results 

We describe here the results of our fits and compare 
them with other investigators' efforts. The fits arc shown 
in Figs. 3-17, and listed and compared with others' in 
Table 3. 

3.6.1. PSR B0301+19 (Fig. 3) 

The two other RVM fits at frequencies similar to or be- 
low ours, NV82 and BCW91, obtain substantially different 
results from our fit. All three yield similar values of the 
ratio sina/sin/?, as expected from Eq. 5. We find that 
the other two fits share similar /3 and appear to have a 
related by tt — a^v ~ olbcw although we believe that we 
have reduced all results to common definitions. Note that 
this transformation would represent a switch from an inner 
to an outer line of sight trajectory. The disparities high- 
light the great difficulty in determining a unique a and (3 
from an RVM fit for these two highly covariant parame- 
ters. The HX97a,b fit is at a much higher frequency and 
not surprisingly has quite different parameters. 

The high quality of our position angle data, representing 
0.3(i of integration at Arecibo Observatory, has enabled us 
to fit farther into the pulse profile wings than could the 
earlier RVM investigators. Apparently this extra longitude 
range permitted a more robust solution. It is quite inter- 
esting to note that our results are much closer to the two 
E/G investigations, which also are similar to each other. 

3.6.2. PSR B0525+21 (Fig. 4) 

The results in Table 3 appear rather scattered at first 
glance. However, We agree with BCQ to within a few a, 
and the HX97a,b 1.41 GHz result would also agree if one 
took the complement of anx (switching from an inner to 
an outer line of sight trajectory). However, the two E/G 
results and the NV82, 0.43 GHz RVM fit place a far away, 
at or above ~ 155°. We favor our result, which is based 
on 1.9h of integration at Arecibo, again permitting a fit 
farther into the profile wings. 

3.6.3. PSR B0656+14 (Fig. 5) 

The low-level polarized emission between (18 — 30)° lon- 
gitude provides part of the critical second half of the classic 
S-shaped RVM curve. It is this extra longitude coverage, 
resulting from a deep 1.8ft. integration, that enables the fit 
to converge to a unique value. Note that our fit moves (/)o, 
the longitude of symmetry, to the end of the main pulse 
component, indicating that the main component and the 
trailing low-level emission constitute the opposite sides of 
a cone. There are no other RVM results on this pulsar. 
Our rather large error bars encompass the E/G results. 

3.6.4. PSR B0823+26 (Figs. 6,7) 

Both R93 and W99 have labelled the main component a 
core, based primarily on its polarization properties. (R93 
called the whole profile T, while W99 labelled it St based 
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also on newer high-frequency observations.) The pulsar 
presents a main pulse plus a postcursor at (j) ^ 40° and 
a very weak interpulse near (j) ~ 180°. An apparent or- 
thogonal mode switch between cj) = (—180 and —7)° and 
(j) — (+140 and +190)° was included in the fit by adding 
90° to the position angle over that longitude range, while 
a region of competing orthogonal modes was unweighted 
in the <j> = (—7 to —1)° range. We present two sets of 
results - the preferred one to the full main / postcursor / 
interpulse component range, and another one to the first 
two components alone (see Table 4). The two results dif- 
fer by several times their formal errors, a point to which 
we shall return below with pulsar B0950+08. Note also 
that our two results place a in (barely) opposite hemi- 
spheres. This seems rather unimportant until one realizes 
that the negative /? of both fits then implies an inner line 
of sight trajectory in one fit, but an outer trajectory in the 
other (see Section 1.1 and Table 2 for further discussion). 
However, the distinction between inner and outer is less 
important for equatorial trajectories such as these. 

Our results comport with the other three RVM fits on 
this pulsar at similar frequencies: NV82, BCW91, and 
HX97a,b, although error bars on their results are gener- 
ally quite large or not given. (A 0.43 GHz BCW91 fit 
with small error bars is also in quite good agreement with 
our 1.4 GHz results.) Both of the empirical/geometrical 
main-pulse results lie close to ours. LM88 also derive E/G 
results for the interpulse, as shown in Table 3. In addi- 
tion, LM88 crudely fitted the RVM to a combination of 
0.43 and 1.4 GHz position angle data for all three pulse 
components, with the result that a = 100°, /3 = —5.2° (no 
uncertainties given). The fourth panel of Figs. 6 and 7, 
which shows the postion angle calculated from many work- 
ers' results atop our data, demonstrates that our fit yields 
a better representation of our data than does any other 
published result. Our careful treatment of uncertainties 
and biases in the extensive low-level emission also favors 
our result. 

It is clear from our fit that the interpulse represents 
emission from the opposite pole having a2 = 81?1, so 
02 ^ C ^ oi2 = 14?8. This large impact parameter helps 
to explain why the interpulse is so much weaker than the 
main pulse. The postcursor is then emitted from the vicin- 
ity of the primary magnetic pole but at higher impact pa- 
rameter than the main pulse. 



3.6.5. PSR 80950^ 



(Figs. 8, 9, 10) 



In our preferred fit (Fig. 8), we use data from all lon- 
gitudes, with the exception of the ranges (j) = (—10 to 
+ 15)°. We also account for the orthogonal mode emission 
between (j) = (25 to 175)°. The resulting fit appears quite 
robust across virtually the whole fitted longitude range. 

Our result is significantly at odds with all the others 
(including both RVM and E/G determinations), which are 
rougly consistent with each other. While our solution sug- 
gests a two-pole orthogonal rotator like B0823+26, the 
other fits (including an LM88 RVM fit to data stitched to- 
gether from several sources yielding {a, (3) = (170,4.6)°) 
place a « 170° with both pulse and interpulse emission 
from a single pole. [Narayan & Vivekanand (1983b) and 
Gil (1983) elaborated a one-pole model for this pulsar.] 

Our result indicates that the second pole has a2 — 74? 6 



and (32 = 52?9. Hankins & Cordes (1981) showed that the 
main pulse-interpulse longitude separation remains con- 
stant between 40 and 4850 MHz, as expected of a two-pole 
emitter. However, the emission between the pulse and in- 
terpulse components is then particularly puzzling, as it is 
emitted even farther from the two magnetic poles than are 
the pulse and interpulse themselves. 

Our two other solutions select data over narrower lon- 
gitude ranges to test the consistency of our fits, (see Ta- 
ble 4 and Figs. 9 and 10). First, we unweight data in the 
vicinity of the main pulse, leaving what appears to be a 
fine fit across the interpulse and elsewhere (Fig. 9). The 
angles a, /3 and (j)o change significantly. The second ad- 
ditional fit unweights all data near the interpulse. Note 
that this "main pulse" fit looks quite reasonable through- 
out the fitted range (see Fig. 10), yet a and (3 move from 
our preferred fit by several times the formal errors. 

While both of these latter fits appear reasonable over 
their fitted ranges, they both fail elsewhere. Only the 
first, full longitude range fit, conforms well to observed 
position angles at all longitude ranges. We believe that 
our high-quality data, coupled with our careful treatment 
of measurement uncertainties in the fits, provide a more 
accurate result. For comparison, we superpose some of the 
other workers' results, along with ours, onto our data in 
the bottom panel of each of Figs. 8, 9, and 10. Note that 
only our full fit (Fig. 8) matches well the overall slope 
of position angle with longitude over the longitude range 
~ —160 to —20°. (All results including ours have some 
trouble matching the left edge of the interpulse, presum- 
ably as a result of orthogonal mode competition.) 

We are confident that our full range fit is the best one. 
However, our own three somewhat inconsistent solutions 
help to illustrate the pitfalls of RVM fitting. Since most 
other pulsars' emission is detected over a much narrower 
longitude range, the resulting solutions could be expected 
to be no more representative of the actual situation at 
the pulsar than are the two restricted range fits discussed 
here. Presumably orthogonal mode mixing and/or other 
systematic effects cause these inconsistencies. Clearly the 
formal uncertainties generally underestimate the errors, as 
was also the case for pulsar B0823+26. 

Similarly, the fact that other RVM investigators con- 
verged on a much different solution than we did (indicat- 
ing a one-pole emitter rather than our two-pole orthog- 
onal rotator) shows how sensitive the fitting process is to 
low-level emission which must be properly weighted. 

3.6.6. PSR B1541+09 (Figs. 11, 12) 

This pulsar has an unusually wide profile, which extends 
over 100° in longitude. Unfortunately a large region, from 
^ — 10 to ^ +35° longitude, must be excluded in order 
to achieve convergence. Note that within the context of 
the RVM model, the only way to explain the reversal of 
sign of (hp/d4> within the profile is to posit an outer line 
of sight trajectory, with the extremum in ijj occuring at 
a longitude that is some distance from (f>Q (Narayan & 
Vivekanand 1982). Indeed our fit, after excluding the 
above noted range, conforms to these conditions. How- 
ever, the full range of observed position angles cannot be 
fitted to any possible set of parameters in this scenario. 
More likely, emission in the excluded range represents the 
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combination of two RVM orthogonal modes with relative 
amplitude varying across the excluded zone. Backer & 
Rankin (1980) show orthogonal mode competition is im- 
portant at lower longitudes than these at 430 MHz; un- 
fortunately there are apparently no 1.4 GHz orthogonal 
mode analyses. 

There are no other published RVM fits, but the two E/G 
results (which are close to one another) are significantly at 
odds with ours. However, the position angles derived from 
both of their results bear no resemblance to the observed 
position angles, which leads us to doubt their efficacy, at 
least at our frequency of 1.4 GHz (see the bottom panels 
of Figs. 11 and 12). 

Proceeding under the assumption that some or all of 
the data at longitudes below ^ 35° must be unweighted, 
we were able to obtain two convergent fits. In our pre- 
ferred fit (see Fig. 11), the longitude range (—8, +38)° is 
unweighted, but lower longitudes fit quite well onto the 
RVM. In the second fit, all longitudes below 30° are un- 
weighted. The two fits yield rather different results for 
a, (3, and 00 . 

In our fits, the inflection point at ^ 52 or ~ 62° longi- 
tude represents the closest approach of the magnetic axis 
and the line of sight. It is interesting that the very strong 
circular polarization (among the strongest observed in any 
pulsar) occurs at longitudes similar to our excluded range, 
suggesting an unusual emission component much like a 
core except that it is not at the symmetry point, as ex- 
pected. W99 reviewed all available multifrcquency profile 
measurements and concluded that this pulsar is a classical 
triple (T), also identifying the highly circularly polarized 
component with the core. The current fit differs from their 
result in that it finds the symmetry center after the highly 
circularly polarized component, rather than in it. The 
leading component is then even farther from the symme- 
try center, rather than lying symmetrically opposite the 
trailing component. 

3.6.7. PSR B1839+09 (Fig. 13) 

Our result with a nearly 90° and (3 near 2° agrees very 
well with the E/G analyses of R93 and LM88, in our fit 
encompassing 80° of longitude. The only other RVM fit, 
from BCW91, did not achieve statistically significant re- 
sults. Note that we unweighted the orthogonal mode com- 
petition starting at longitude 4°. 

3.6.8. PSR B1915+13 (Fig. 14) 

For this pulsar, we fit to all data within 60° of longitude 
centered around the main pulse. The classical RVM sweep 
appears to fit the data well. Comparing our results to oth- 
ers, the BCW91 error bars on a and f3 barely overlap with 
ours. As our fit is based on significantly more data, it is 
probably more reliable. The agreement with HX97a,b at 
4.85 GHz is poorer, which is not surprising given the large 
frequency difference. The E/G results of R93a,b are the 
closest solutions to ours. 

3.6.9. PSR B1916+14 (Fig. 15) 

We fit the central 13° of longitude. The only other RVM 
fit, from BCW91, has very large uncertainties which en- 
compass our fit. Of the two E/G results, LM88 find a and 



P near ours but the RVM curve derived from their parame- 
ters deviates significantly from our data, while R93a,b ex- 
hibits the opposite situation: her a and /3 are far from ours 
but the RVM curve derived from her parameters matches 
our data well. The latter situation results from the strong 
covariance of the two parameters. 

3.6.10. PSR B1929+10 (Figs. 16, 17) 

For our preferred fit (Fig. 16), we fit all 360° of longi- 
tude, except that we unweight the main pulse region from 
— 17° to -1-10°. The excluded region could not be fitted sat- 
isfactorily with the RVM curve, especially when trying to 
fit the interpulse and other off-pulse emission simultane- 
ously with the main pulse. The single-pulse polarimetry of 
Rankin & Rathnasree (1997) [RR97] shows that indeed a 
competition between two orthogonal modes leads to posi- 
tion angle distortions in this region. Because the pulsar is 
so strong at these longitudes, our weighting scheme would 
otherwise cause this region to dominate the fit. We also 
placed a 90° orthogonal mode offset from longitude —60° 
to —15°, which is also evident in the single-pulse displays 
of RR97. Our second fit eliminates the main pulse entirely, 
(see Table 4 and Fig. 17)), and yields results roughly sim- 
ilar to the first. 

Our a of approximately 35° and (3 ^ 26° seems to be 
another shot at what seems a scatter of results for this pul- 
sar. Other results not listed in Table 3 include the follow- 
ing RVM fits: LM88, using the 430 MHz data of Rankin & 
Benson (1981), found (a, /3) = (15, 7.5)°, Phillips (1990) 
[P90] measured (a,/3) = (35 ± 4,21 ± 3)° at 1665 MHz 
and (a,/3) = (30 ± 2, 20 ± 2)° at 430 MHz, and RR97 
found {a, (3) = (31,20)° at 430 MHz. Note that most of 
the RVM fits, including ours, yield roughly similar results, 
which indicate that the pulse and interpulse are emitted 
from nearly opposite sides of a wide, hollow cone near the 
+Q. spin axis. While one can find the colatitude and im- 
pact parameter of the other magnetic pole as a2 = 144? 
and P2 = — 82?5, these quantities are not important in a 
one-pole emission model. The significantly nonzero value 
of 00 indicates that the brightest emission does not occur 
at closest approach of the line of sight and the magnetic 
axis. Aside from his choice of emission centered on the 
— f2 spin axis, we are in agreement with the rough sketch 
of the emission geometry shown in P90. 

R93, using E/G arguments, also found [a, (3) = 
(18, 11.6)° ff B1929+10 is interpreted as a cT-type pulsar, 
but suggested that it might be an orthogonal rotator if it 
is of the T class. While the RR97 RVM fit discussed above 
also suggests a moderately tipped dipole, the authors ar- 
gue that RVM fits lead to spurious results for this pulsar, 
and they favor an orthogonal rotator {a, (3) ^ (90,0)°, 
dual pole emission model, again on the basis of E/G con- 
siderations. We still prefer our RVM results, as we believe 
it is unlikely that such a beautiful fit obtained over such 
a wide longitude range could be so seriously corrupted by 
systematic effects as to render it invalid. 

4. CONCLUSIONS 

We have fitted pulsar polarization data to the rotating 
vector model (RVM) in order to find geometrical parame- 
ters. We were careful to account properly for the bias in 
L at low S/N and for non-Gaussian statistics of position 
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angle data, so that the low S/N data could be appropri- 
ately weighted in the fits. This enabled us to correctly 
include low-level emission arising from weak pulse compo- 
nents as well as the wings of main components in the fits. 
We are not aware of any other RVM fits that include such 
considerations. 

We succeeded in finding convergent RVM fits on ten of 
the pulsars observed by Weisberg et al. (1999). In all 
such cases, we always converged to our same solution re- 
gardless of our choice of initial parameters gleaned from 
previously published results. Yet we were surprised that 
our fitting procedure converged for relatively few of the 
many pulsars for which we had high-quality data. Pre- 
sumably this was due to competition between orthogonal 
emission modes or other systematic effects. Because of the 
high covariance between the two important fitted parame- 
ters a and /3, even slight deviations of position angle data 
from the nominal RVM model could prevent convergence 
to a unique solution. Wide longitude ranges of emission 
were generally required in order to achieve unique, con- 
vergent fits. As a result, pulsars with interpulses and/or 
other off pulse emission are favored targets for RVM stud- 
ies. Nevertheless, we show that in two such cases (pulsars 



B0823-f26 and B0950+08), our fits to different longitude 
ranges yield results differing by several times their formal 
errors. Consequently it is clear that systematic effects are 
present at a significant level. 

We created transformation equations to convert other 
published RVM and empirical-geometrical (E/G) results 
to a common definition, in order to compare them with our 
results and to facilitate future comparisons. For some pul- 
sars, all workers' results are similar while for others there 
is a wide range of solutions. No clear pattern emerges from 
these comparisons. 

We are currently in the process of separating orthogo- 
nal modes in single pulse polarization data (Weisberg & 
Everett 2001), in order to to determine if RVM fits to 
separated emission modes yield results different than the 
average pulse results presented here. 
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Fig. 1. — In the rotating vector model (RVM, see Eq. 1), the magnetic pole colatitude a, line of sight colatitude (J, and line of sight impact 
parameter /3 are defined in terms of the positive angular velocity vector +Q and the observable magnetic pole fi. These quantities are shown 
here in the plane containing Q and fi. Note that these definitions are the same for all possible a; i.e., for (0 < a < tt) radians. 
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Fig. 2. — X surface as a function of a and /3 for pulsar B0301+19 Note that there is a long, deep, and virtually flat— bottomed chasm 
defined by the observed value of position angle sweep rate dip/dcj)^ via Eq. 5. 
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Longitude, [deg] 

Fig. 3. — Pulsar B0301+19: (a): Total (/), linearly polarized (L), and circularly polarized {V, dashcd-dotted) , flux densities; (b): Observed 
and fitted position angle of linear polarization i/i with error bars representing measurement uncertainties; and (c): Position angle residuals, 
normalized by a^ . The fit was allowed to extend over the full displayed longitude range. 
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Fig. 4.— Pulsar B0525+21: Sec Fig. 3 caption for details. 
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Fig. 5. — Pulsar B0656+14 Sec Fig. 3 caption for details. 
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Fig. 6. — Pulsar B0823+26. Our preferred full longitude fit. For an explanation of (a) — (c), sec Fig. 3 caption. Note that a 90° orthogonal 
mode switch is inserted in the longitude range (—180, —7)° and (+140, +180)°, and that longitudes between (—7, —1)° are unweighted in the 
fit. Panel (d) shows the RVM curves from our and other investigators' work. Orthogonal mode jumps included in our fit and shown in (b) and 
implicitly in residuals (c) are not included in (d), in order to simplify comparison with other authors' work. Also, since other investigators 
generally did not provide their fitted position angle and longitude offsets tpo and ij>0 > we use our values of these two parameters in calculations 
of their RVM curves. At the resolution of panel (d), the NV82 RVM is indistiguishable from ours, so we do not plot it. 
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Fig. 7. — Pulsar B0823+26. Main pulse through postcursor fit. For general explanation of this figure, see Fig. 3 and Fig. 6 captions. Note 
that a 90° orthogonal mode switch is inserted in the longitude range (—20, —7)° and that longitudes between (—7, —1)° are unweighted in 
the fits shown in (b) and (c). At the resolution of panel (d), both the BCW91 RVM and the NV82 RVM are indistiguishable from ours, so 
we do not plot them. Also, the HX97 RVM curve overlays the LM88 E/G curve, so we do not plot the HX97 RVM fit. 
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Fig. 8. — Pulsar B0950+08. Our preferred full longitude fit. For general explanation of this figure, see Fig. 3 and Fig. 6 captions. Note 
that a 90° orthogonal mode switch is inserted in the longitude range (25, 175)° and that longitudes between (—10, +15)° are unweighted in 
the fits shown in (b) and (c). At the resolution of panel (d), the R93 E/G curve is indistiguishable from the LM88 E/G curve, so we do not 
plot the R93 E/G curve. 
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Fig. 9. — Pulsar B0950+08. The interpulse is weighted in this fit, with the main pulse ((j> = (—90, 90)°) unweighted. For general explanation 
of this figure, see Fig. 3 and Fig. 6 captions. Note that a 90° orthogonal mode switch is inserted in the longitude range (25, 170)° and that 
longitudes between (—90,90)° are unweighted in the fits shown in (b) and (c). We omit the same curves as in Fig. 8. 
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Fig. 10. — Pulsar B0950+08. Main pulse weighted, with the intcrpulse unweighted. For general explanation of this figure, see Fig. 3 
and Fig. 6 captions. Note that a 90° orthogonal mode switch is inserted in the longitude range (25,170)° and that longitudes between 
(—180, —90)°, (—10, 15)°, and (90, 180)° are unweighted in the fit shown in (b) and (c). We omit the same curves as in Fig. 8. 
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Fig. 11. — Pulsar B1541+09. Our preferred wide longitude fit. For general explanation of this figure, see Fig. 3 and Fig. 6 captions. 
Longitudes between (—8, +38)° are unweighted in the fit shown in (b) and (c). The R93 curve in panel (d), which appears to have the wrong 
slope, corresponds to {a, (3) = (175, 0), which represents a line of sight traverse right across the magnetic pole. Even a slight nonzero (3 would 
add a region of inverted slope near (f> = (f>o, which would better match our data. 
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Fig. 12. — Pulsar B1541+09. Narrower longitude fit. For general explanation of this figure, see Fig. 3 and Fig. 6 captions. All longitudes 
below 30° are unweighted in in this alternative fit shown in (b) and (c). See Fig. 11 for a discussion of the R93 curve in (d). 
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Fig. 13. — Pulsar B1839+09. Sec Fig. 3 caption for details. Longitudes between (+4, +10)° are unweighted in the fit. 
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Fig. 14. — Pulsar B1915+13. Sec Fig. 3 caption for details. 
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Fig. 15. — Pulsar B1916+14. For a general explanation of this figure, see Fig. 3 and Fig. 6 captions. At the resolution of panel (d), the 
R93 E/G curves is indistiguishable from the BCW91 curve, so we do not plot the R93 E/G curve. 
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Fig. 16. — Pulsar B1929+10. Our preferred full longitude fit. For general explanation of this figure, see Fig. 3 and Fig. 6 captions. Note 
that a 90° orthogonal mode switch is inserted in the longitude range (—60, —15)° and that longitudes between (—17, +10)° are unweighted 
in the fits shown in (b) and (c). Our preferred RVM curve and most other published RVM fits provide the best match to the data. At the 
resolution of panel (d), the NV82 RVM, RR97 RVM (their fit), and BCW91 RVM curves are indistiguishable from our curve, so we do not 
plot them. Also, the LM88 RVM (their fit) curve is indistinguishable from the LM88 E/G (their empirical calculation) curve, so we do not 
plot the LM88 RVM curve. 
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Fig. 17. — Pulsar B1929+10. Entire main pulse region unweighted. For general explanation of this figure, see Fig. 3 and Fig. 6 
captions. Note that a 90° orthogonal mode switch is inserted in the longitude range (—60,-15)° and that longitudes between (—17, +10)° 
arc unweighted in the fits shown in (b) and (c). Our RVM curve and most other published RVM fits provide the best match to the data. We 
omit the same curves as in Fig. 16. 
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Table 1 
Dictionary and Conversion Table from Earlier Work for Geometrical Beam Parameters. (See Eqs. 1 and 2 

AND ensuing discussion.) 



Investigators 
and 

Method 


Convention 
Problem? 


Colatitudc of Observable Magnetic Pole, a. 

w.r.t. (+n) Spin Axis 


Impact Parameter of Line of Sight 
w.r.t. Observable Magnetic Pole 


Conlined 

to First 

Quadrant? 


Relation to 
our a 

[c^OTiqii^al -^ our Cx) 


Symbol in 

Original 

Paper 


Hclation to 

our /3 
(Symbolariqinal ~* our /3) 


Current Work 
(RVM) 


no 


no 


a — > Q 


f:i 


0-^0 


NV82 
(RVM) 


no 


yes 


If {dip/dcp and Pnv) 
both have same sign, then 

TT — QiVV — * Ct- 

Otherwise, c^nv — ^ «• 


f^NV 


If dip/d(t> > 0, then 

-|/3jvvl -» /3. 

If dif>/d4, < 0, then 

\Pnv\ -^/3. 


LM88 
(E/G) 


n/a 


yes 


As the sign of f^^M was not 
published, either aLM — * ct 

or TT — aLM — * CK are possible. 

We selected one based on con- 
sistency with our fits. 


f^LM 

(magnitude only 
was published) 


If dip/d(t> > 0, then 

If dip/d<p < 0, then 
|/3lm| -^13- 


B0W91 

(RVM) 


yes 


no 


TT — ascw — * a- 


<^BCW 


<^BCW -^ —1^ 


R90,93a,b 

(E/G) 


n/a 


yes 


If {d'iP/d4) and 0h) 

both have same sign, then 

TT — an — » a. 

Otherwise, an — > a. 


fiB. 


If dip/d(t> > 0, then 

-|/3b| ^ /3. 

If di/j/d0 < 0, then 

I/3h|^/3. 


HX97a,b 

(RVM) 


yes 


no 


TT — \aHx\ -^ a. 


ohx 


it a/fx > 0, then 

(7HX — * — /3. 

Otherwise, cr/fx -^ /3-^ 



Note. — (a) In one case, (PSR B0950+08 at 1.41 GHz), the resulting (^ would then be greater than tt. Consequently, 2\aHx\ — <^HX -^ /3 
must instead be used. 



Table 2 
Inner and Outer Line of Sight Tra,iectories 



Sign of 
Slope" 

dip 1 

d<p \^°-^ 


Sign of 
Slope'' 

dip' \ 
d<P 1™°^ 


Sign of Impact Param. 

of Line of Sight w.r.t. 

Obs. Mag. Pole, /? 


Colatitude of 
Observable Mag- 
netic Pole, a 


Line of Sight Trajectory 

w.r.t. 

Observable Magnetic Pole 


Negative 


Positive 


Positive 


<7r/2 
>7r/2 


Outer (Equatorward)^ 

Inner (Nearest [— f7] Spin Poleward) 


Positive 


Negative 


Negative 


<7r/2 
>7r/2 


Inner (Nearest [+0] Spin Poleward) 
Outer (Equatorward)^ 



Note. — (a) The quantity -rrlmax is the position angle sweep rate in the observers' convention (increasing counterclockwise on the sky), 

while (b) -ji-\max is the sweep rate in the (opposite) RVM convention, (c) "Equatorward" indicates that the line of sight is opposite the spin 
pole lying nearest to the observable magnetic pole. 



Table 3 
Geometrical angles a and (i from our work and other investigators.' All angles have been converted to a 

COMMON DEFINITION USING TABLE 1. 



Is2 

00 



Pulsar 


R93a,b 


O 


nr adopted results 


NV82 


LM88 


BCW91 




R93 


a,b 


HX97a,b 






Class. 




(1.42 GHz) 




(0.43 GHz) 


(0.4 GHz) 








1 GHz 














(RVM) 




(RVM) 


(E/G) 


(RVM) 




(E/G) 


(RVM) 








a 


/3 


4>o 


xi 


a (3 


a 


/3 


V (GHz) 


a 


/3 


a 


/3 


V (GHz) 


a 


/3 


0301+19 


D 


162.4 
±11.8 


0.96 
±0.63 


-3.71 
±0.006 


15.6 


110 3 

±11 


148.1 


1.8 


1.42 
0.43 


69 
±16 

61 
±30 


2.9 
±0.3 

2.8 
±0.4 


150 


1.7 


4.85 


85.7 
±15.0 


13.9 
±15.0 


0525+21 


D 


116.8 
±4.6 


-1.50 
±0.08 


-5.84 
±0.02 


7.9 


> 160 -0.6 


156.8 


-0.7 


1.42 
0.43 


134 

±10 

140 

±9 


-1.2 
±0.2 

-1.0 
±0.2 


159 


-0.6 


1.41 
1.71 
4.85 


44.2 

±60 
113.7 

±80 
162.2 

±70 


-1.3 
±0.2 

-2.0 
±1.0 

-0.6 
±2.0 


0656+14 


rpa 


29 
±23 


8.9 
±6.1 


14.9 
±0.7 


2.6 




8.2 


8.2 








30 










0823+26 


St 


98.9 
±0.7 


-3.03 
±0.01 


1.33 
±0.01 


8.6 


98 -4 


76.9 


-1.1 


1.42 


89 
±600 


-3 

±1 


84 - 


-1.9" 


1.41^= 


90.2 
±25 


-1.0 
±0.2 


(Main / 


















0.43 


101 
±1 


-3.2 
±0.3 






1.71= 


97.7 
±80 


-1.7 
±1.5 


Postcursor) 




























4.85^= 


82.5 
± 15 


12.1 
± 15 


Interpulse 














±26.1 


5.6 


















0950+08 


Sd 


105.4 
±0.5 


22.1 
±0.1 


-11.9 
±0.2 


2.2 


170 5 


174.1 


4.2 


1.42 


174 
±30 


2.5 

±15 


168 


8.5 


1.41 


179.3 

±15 


0.4 

±15 


(Full) 


















0.43 


174 
±90 


2.5 
±40 






4.85 


153.1 
±90 


4.1 
±5 


1541+09 


T 


131.0 

±5.67 


-20.23 
±2.26 


62.5 
±0.37 


1.5 




174.2 


-0.8 








175 


0.0 








1839+09 


T 


86.1 
±11.4 


2.3 
±0.04 


0.34 
±0.02 


10.3 




90 


2.9 


1.42 


146 
±400 


1 

±8 


97 


1.4'' 








1915+13 


St 


73 

±19.4 


5.4 
±0.5 


6.8 
±0.1 


4.2 








1.42 


86 

±24 


5.4 
±0.2 


68 


6.6 


4.85 


89.3 

±12.5 


9.3 

±15 


1916+14 


'!■(■!) 


118.0 
±33.4 


-1.0 
±0.33 


0.38 
±0.01 


8.9 




124.5 


-2.1 


1.42 


167 
±100 


-0.3 

±0.2 


101 


-1.3 








1929 + 10" 


T(?)/ 
cT 


35.97 
±0.95 


25.55 
±0.87 


-19.62 
±0.2 


2.45 


35 23 


6 


4 


1.42 

1.42<= 

0.43 


27 

±4 

150 

±10 

25 
±2 


16 

±3 
-3 

±2 
16 

±2 


90^ 

18/ 


41.8 
11.6 








Interpulse 
























88 











< 



^ 



0) 



Note. — (a) W99 suggest St; (b) Here we choose opposite sign as R93a,b, as required by Eq. (5); (c) Main pulse only; (d) See text for additional 
RVM fits by LM88,P90, and RR97; (e) Alternate BCW91 fit; (f) First Rankin fit is under assumption of T classification; second is for cT. 
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Table 4 
Multiple convergent fits to pulsars with multiple components. 



Pulsar Components fitted Longitude 

limits" 



Our results 

/3 



Po Xu 



0823+26 



full 



±180 98.9 -3.03 1.33 8.6 

± 0.7 ± 0.01 ± 0.01 



main / postcursor -20 to 45 85.8 -3.08 1.21 8.6 

± 3.2 ± 0.01 ± 0.02 



0950+08 



full 



±180 105.4 22.1 -11.9 2.2 

± 0.5 ± 0.1 ± 0.2 



main pulse -90 to +90 126.2 19.6 -10.8 1.3 

± 4.3 ±1.0 ± 0.3 



interpulse 



-90 to -90 



109.6 13.3 
± 1.4 ± 0.9 



-37.6 2.3 
±5.3 



1929+10 



full 



±180 35.97 25.55 -19.62 2.45 

± 0.95 ± 0.87 ± 0.2 



interpulse -180 to -40 34.3 28.8 -20.2 2.5 

35 to 180 ± 0.86 ±0.23 ± 1.2 



Note. — (a) See text and figure captions to find ranges of unweighted 
longitudes within these limits. 



